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Abstract 

We study the effects of cosmic strings on structure formation in open uni- 
verses. We calculate the power spectrum of density perturbations for two 
class of models: one in which all the dark matter is non baryonic (CDM) and 
one in which it is all baryonic (BDM). Our results are compared to the 1 in 6 
IRAS QDOT power spectrum. The best candidates are then used to estimate 
fj,, the energy per unit length of the string network. Some comments are made 

on mechanisms by which structures are formed in the two theories. 
PACS Numbers : 98.80.Cq, 95.35+d 
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The cosmic string scenario has generated much interest over the past few years as a 
promising candidate for structure formation. Following a primordial phase transition, line- 
like concentrations of energy could form in certain Grand Unified Theories [1]. These would 
then evolve into a scaling regime, in which all correlation functions should be expressed 
in terms of dimensionless ratios of length scales over the horizon. The long-string energy 
density would evolve as 



perturbing the background radiation and matter density, leading to the formation of astro- 
physical objects and anisotropics in the microwave background. 

Albrecht & Stebbins have performed an analysis of string-seeded perturbations in non- 
baryonic Cold Dark Matter (CDM) and Hot Dark Matter (HDM) [2]. They have developed 
a model (form here on the AS model) which not only allows for different possibilities for 
string network evolution and matter content of the universe but also includes the effect of 
compensation. They have found that, in a CDM universe, strings produce too much power 
on small scales relative to large scales, more so than in an adiabatic CDM cosmology. They 
propose however that in a HDM cosmology the freestreaming by the neutrinos will partially 
suppress perturbations on small scales to a level which might be consistent with the data. 
Recent analysis of large scale surveys have indicated that the spectrum of mass fluctuations 
tends to flatten out on scales larger than in standard (0 = 1) CDM or HDM models, and 
consequently the effect of HDM on string perturbations would not be effective on large 
enough scales. Thus strings and HDM seems to have some of the same problems as its CDM 
counterpart. 

Although there are strong theoretical reasons to favour an = 1 universe, as yet there is 
no strong observational evidence for it (see [4]). On the other hand, the extension of standard 
cosmological scenarios to open universes has proven to be quite promising; CDM adiabatic 
and topological defect models do not have excess small scale power if one resorts to low fi/j 
universes. This renewed interest in open universes has led us to reanalize the AS model. Due 
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to the large number of choices we can make (string network model, type and density of dark 
matter, H 0} ionization history, etc) we have opted to present the results of what we think 
are representative models. We will consider standard recombination, H = 50Kms _1 Mpc _1 , 
universes where all the dark matter is either non-baryonic or baryonic (BDM) and we have 
chosen one model for the network evolution. In a future publication we shall present a more 
detailed description of parameter space. 

We shall briefly illustrate the formalism we are using. One is interested in calculating the 
power spectrum of density perturbations, defined as (2ir) 3 P (k)8^ (k — k') = (6(k)6(— k')). 
The AS model uses the form 

P(k) ~ I67r 2 (G7i) 2 / |G M (MoY)| 2 .W/a)<V (2) 

where: 

1) The "structure function": as proposed in [2], 

is the time integral of the two-point correlation function of the u p + 3P" part of the long 
string stress energy tensor. We must emphazise that it is specific to cosmic string physics; 
it has a k~ 2 behaviour on scales smaller than the long string coherence scale, characteristic 
of sheets, and a k° behaviour on large scales, typical of uncorrelated point-like objects. The 
different parameters are £ = (/^//i)^, the long string coherence or "curvature" scale x, the 
r.m.s. velocity of strings f3 and the amount of small scale structure E. 

2) The string network: we are assuming that the main source of perturbations are the 
long strings [5] . There are no numerical simulations of string networks in a curvature domi- 
nated universe and all current efforts to find reliable analytical models have focused on flat 
universes. We have modeled the string network in three different ways. Having determined 
the differences to be unimportant we have chosen one model for subsequent calculations. 
Model 1 is the original idealized model which has the correlation length scaling with the 
horizon, the X model in [2]. Model 2 is a good fit to the flat-spacetime simulations and is 
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given in [6]; we extrapolate to an open universe using the parameters calculated deep in the 
radiation era in [7] [8]. In this case, when curvature starts to dominate, strings will slow 
down and one finds that the network will deviate from scaling. We find the same sort of be- 
haviour in model 3 proposed in [3], motivated by large-N non-linear sigma model dynamics 
and by nematic liquid crystals where evolution is friction dominated. Here £ ~ 772 in the 
free-expansion phase, so one can interpolate between the scaling and nonscaling behaviour 
with x = One finds that, on the scales we are interested, deviation from scaling 

will have no significant effect. In fact all these models give very similar powerspectra albeit 
with different normalizations. We have chosen to work with model 2 so x = 2(, E/3 2 = 1.2, 
and k c = Att/t]. 

3) Green's function: Gm is the Green's function of the matter density perturbation which 
undergoes gravitational collapse at late time. The equations are solved in the synchronous 
gauge (all perturbations are seeded inside the horizon so we don't have to deal with gauge 
issues) and k is the eigenvalue of the Laplacian on a hyperbolic space. As in AS we are 
correcting for compensation by including a factor of 1/(1 + (k c /k) 2 ) which gives the correct 
energy momentum behaviour on large scales. 

The evolution equations for the Green's functions are integrated numerically. In the CDM 
case we have checked that they match, for = 1, the analytic fit of [9] and, for variable 
and h and no sources, the pure CDM transfer function of [10]. For BDM, following [11], 
we have solved the fluid equations in the tightly coupled regime until this approximation 
breaks down. We then calculate explicitly the effect due to damping during recombination 
and finally match onto post-recombination solutions. The "Sachs- Wolfe" term is simplified 
by discarding the traceless part of the metric, and we model the evolution of the (standard 
recombination) ionization fraction accurately by integrating the equations derived in [12]. 

We can fit our results with the function: 





where the a's can be read off from the table (k in units of h 2 Mpc 1 ). 
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a 2 


a 3 


ci4 


a 5 


bGfie 


&BDM = 


.2 


143.7 


23 


21 


17.5 


.00031 


9.6 


&BDM = 


.4 


190.0 


11.5 


11 


12.55 


.0004 


4.6 


&CDM = 


.1 


195.2 


47 


44 


15.5 


.0002 


9.9 


&CDM = 


.3 


168.1 


15.7 


14.7 


5.17 


.00037 


4.3 



TABLE I. Parameters of four models normalized, to fit the 1 in 6 QDOT IRAS power spectrum; 



b is the bias factor and GfJ, is in units of 10 
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FIG. 1. A plot of the mass fluctuation power spectrum for 4 different models: solid line is 
strings and BDM (0 = .2), dotted line is strings and CDM (0 = .2), long dashed line is strings 
and CDM (0 = 1) and dot-dashed line is adiabatic CDM (0 = 1). They are normalized to agree 
on large scales. 
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As expected, by lowering one can shift the turnover in the power spectra to large 
enough scales. In fact for the CDM models we confirm [4]: one can obtain the lower 
models by simply rescaling, k — > &/fi, a consequence of how late time evolution of the 
network is unimportant on the scales we are interested. In Fig.l we compare two models 
with a standard CDM model and an = 1 strings and CDM model. Note that, unlike 
in adiabatic BDM models, in the strings and BDM model there is no remnant of acoustic 
oscillations, seeing as one is incoherently summing over perturbations seeded at different 
times. 

Although it would be premature to use large scale surveys to make definitive statements 
about our models we have found it instructive to compare four models with the power 
spectrum derived in [13] from the l-in-6 QDOT IRAS redshift survey (Figs. 2 and 3). Clearly 
an ftcDM = -3, h = .5 CDM model continues to have too much power on smaller scales 
(k > 0.07) scales compared to large scales (we have chosen the model in which the turnover 
best agrees with the data), while in an Qcdm = -1 model, although having less power on 
small scales, has a turnover on too large a scale. In a BDM model, the power spectrum 
agrees well with the data for k > .04, but we find that the model that best fits the data on 
these scales (Qbdm = -2) has a turnover on too large a scale. A better fit in this sense is then 
&bdm = -4 [14]. In Table 1. we present an estimate of Gfi which comes from normalizing 
the power spectrum to these scales. Seeing as we are working in linear theory we have opted 
not to normalize to <r 8 where gravitational collapse has most likely already begun. 

Within this model we are unable to make any precise predictions about the cosmic 
microwave background; flat spacetime simulations indicate that all three modes (scalar, 
vector and tensor) contribute on equal footing to the final anisotropy a fact which makes it 
necessary to have a much deeper understanding of the network and its correlation functions. 
In [4] it was pointed out that low multipole moments of the temperature correlation function 
are suppressed; for a CDM, = .2 defect-seeded universe, the spectrum for / < 5 is lower 
than for a flat universe. One would expect, however, that for higher multipoles the spectrum 
is roughly the same in both open and flat models. To reconcile the values of Gfi presented 
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in the table with the one presented in [15], one would need values of b ~ 2.5 — 5. 
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FIG. 2. The power spectrum of two CDM models (O = .1 solid line and fi 
plotted against the 1 in 6 QDOT IRAS power spectrum. 
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FIG. 3. The power spectrum of two BDM models (0 = .2 solid line and i7 



.4 dashed line) 



plotted against the 1 in 6 QDOT IRAS power spectrum. 



The evolution of small wavelength perturbations in BDM is qualitatively different from 



that of perturbations in CDM; collisional (or Silk) damping at around the time of recom- 
bination will suppress baryonic perturbations generated before then (unlike in CDM where 
there is no such effect). However the constant sourcing by the string network will ensure 
that perturbations are generated on small scales at late times (primordial isocurvature per- 
turbations, i.e. perturbations in the radiation field generated at the phase transition, are on 
sufficiently small scales to be ignored). Unlike in HDM, the damping becomes irrelevant in 
a finite time so the tail reflects the k~ 2 behaviour of the structure function on small scales. 

This essential difference between the two models may lead to different features in the 
matter distribution. On small scales in CDM, wakes are competing with an r.m.s. 8m jm 
which has been growing for a long time, beginning in the radiation era. The result is that 
they are not so likely to stand out in a final realization of the matter distribution. In 
BDM however, small scale perturbations have been wiped out, and therefore it is post- 
recombination wakes that are responsible for perturbations on small scales. A possible 
picture is that small scale structures will form in such away as to be organized along large 
scale coherent structures, a realization of the the "radical biasing" suggested in [2]. 

This has bearing on another phenomena of some importance; it has been argued that it 
is likely that a defect-seeded cosmology will be reionized due to early formation of non-linear 
objects, intimately related to the non-gaussian nature of the perturbations. This has been 
verified for texture cosmogonies (where approximately 1% of the matter had formed bound 
objects before a redshift of 50). It is not so clear that one can apply this sort of reasoning to 
string seeded perturbations; the string network tends to be much denser than other defects 
and analysis of both CMBR maps [15] and wake prominence in CDM models indicates 
that string perturbations tend to be much more "gaussian". One expects therefore that 
reionization is less likely in these (CDM) models [16]. There are additional complications in 
the BDM models where on the one hand structure formation can be delayed if the medium 
is ionized but on the other hand wakes can stand out on small scales leading possibly to 
early star formation. We intend to investigate the consequences of late time recombination 
in string models in some detail. 
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In summary, we have looked at string seeded perturbations in open universes and selected 
a few models that are consistent with current large scale surveys; we have considered models 
with either baryonic or non-baryonic dark matter and have tried to identify some of the 
differences in how structures would form on small scales. 
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